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I. INTRODUCTION 

Damping rates of hard particles in a relativistic plasma at high temperature attracted much interest recently p|~^9| . 
For instance, the damping rate of a quark in a QCD plasma is related to imp ortant quantities of the quark-gluon 
plasma, such as mean free paths, thermalization times fl8| , po|| , viscosity |l8| , [2lf , and stopping power P,ffl,g2|-g5| . 
Damping rates also represent simple tests for methods of the finite temperature field theory, probing e.g. the infrared 
behavior jf(| and gauge dependence (l^Jll|,|l^,|l^,|6| ^] of the methods. 

Hard quarks and gluons, i.e. partons with energies E^T, as well as hard electrons and muons in a relativistic QED 
or QCD plasma are damped to lowest order by elastic scattering off the thermal particles mediated by the exchange 
of a gauge boson. Owing to the long range interaction (massless gauge boson) those processes exhibit a quadratic 
infrared singularity. The infrared behavior can be improved by resumming the gauge boson self energy in the "Hard 
Thermal Loop" approximation, leading to an effective gauge boson propagator containing screening effects (Debye 
screening). This prescription is a special case of the Braaten-Pisarski resummation technique p9| ] which produces 
consistent results, i.e. gauge independent results for observables that are complete to leading order in the coupling 
constant. The use of the effective gauge boson propagator, containing the square of the coupling constant g 2 in 
the denominator, reduces the order of the damping rates, naively expected to be proportional to <? 4 , to g 2 , which 
was called anomalous (large) damping At the same time the quadratic infrared singularity is replaced by a 

logarithmic one, reflecting the absence of static magnetic screening in the effective gauge boson propagator. (This 
result led to speculations that these rates might become finite by higher order effects such as a non-perturbative 
magnetic screening mass or self consistent damping [M or that these rates cannot be regarded as direct observables 
but should be modified by transport factors 0-l,Eq].) 

Here we consider the damping of a high energy photon with four momentum P = (po,p) and p = |p| 3> T in a 
relativistic plasma of electrons and positrons or quarks at a temperature much higher than the bare masses of those 
particles. In contrast to partons or leptons the leading order damping mechanism is not caused by elastic scattering 
since photons show no self interaction. To lowest order the photons are damped by Compton scattering and pair 
creation processes shown in Fig.l. Neglecting the fermion masses the corresponding rate suffers from a logarithmic 
infrared singularity. However, this divergence is screened by medium effects again leading to an effective temperature 
dependent fermion mass of the order of gT which cuts off the singularity much more effectively than the bare mass. 
Since the divergence is only logarithmic utilizing a bare propagator it disappears completely taking the effective 
one into account, analogously to the energy loss of an energetic parton in the quark-gluon plasma |7 22 240. Thus 



the introduction of a magnetic screening mass, self consistently damped fermion propagators, or transport factors is 
unnecessary in the case of photon damping. 

Similarly to the energy loss the photon damping rate can be derived consistently by applying the Braaten- Yuan 
prescription |3(]] . For this purpose we decompose the rate into a soft contribution, which requires an effective fermion 
propagator, and a hard contribution, where bare propagators are sufficient. The decomposition is accomplished by 
introducing a separation scale A in the momentum of the exchanged fermion which drops out at the end after adding 
soft and hard contributions if A is restricted by gT « A « T, which is possible in the weak coupling limit j C 1. 
(In the case of partons or leptons the hard contribution can be neglected since the corresponding damping rates are 
dominated by small momentum transfers.) 

The damping rate is defined as the imaginary part of the photon dispersion relation, 7 = —ImE(p), where the 
dispersion relation for a real, i.e. transverse, photon follows from 

pI-p 2 -n T (po,p) = o, (1) 
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where the transverse part of the photon self energy is given by 

1 ( x PiP 3 



Ht (po , p) = 5 ( 5 « - ^ ) n« (po ,P). (2) 
In the case of no overdamping (p <C 7) we get from (Q) 

7= -^-Im n r (p = P,p)- (3) 
2p 

The damping rate is related to the decay rate 1^ by 7 = 1^/2, where T ( i follows from a momentum integration over 
matrix elements and distribution functions (see (||) and @) |pl[ . We will calculate the soft contribution contribution 
starting from ([3]), whereas the hard contribution starting from the decay rate. 

The soft contribution (feo, k < A) follows from the photon self energy shown in Fig. 2, where the blob denotes the 
effective fermion propagator containing the one-loop fcrmion self energy in the Hard Thermal Loop approximation 
p9| . (According to the rules of the Braaten-Pisarski method we do not need an effective vertex or two effective fermion 
propagators at the same time since the external photon is hard.) Cutting through the blob and a bare propagator we 
see that these diagrams correspond to Compton scattering and pair creation via the exchange of a dressed (collective) 
fermion. 

First we will consider the damping of a photon in a hot relativistic QED plasma. The photon self energy shown in 
Fig. 2 is given by 



n„, ' / ^*rh„.S'.g.):,.S*(/0]. (4) 



where S*(K) indicates the effective fermion propagator, the factor of 2 comes from adding the both diagrams in Fig. 2, 
and Q = P — K . It is convenient to use the helicity representation of the fcrmion propagator given in [^,^3|. (The 
effective electron propagator differs from the quark propagator only by the effective mass m 2 = e 2 T 2 /8 compared to 
m q = g 2 T 2 /6.) 

Evaluating the trace over the gamma matrices we get for the transverse part of the polarization tensor 



II T = 2 i e 2 



1 fl-V l + V\ 1 fl + V 1-V 



D+(K)\d+(Q) d-(Q)J D_(K)\d+{Q) d_(Q) 



(5) 



where the functions D±(K) and d±(Q) are given for example by Kapusta et al. |3^] and V — (p • k)(p • <i)/(p 2 qk). 

Using the imaginary time formalism the integration over fco m (r) is replaced by a discrete sum, which can be 
evaluated most easily by adopting the spectral representation for the fermion propagators [ ]32"| , |33| . Using furthermore 
p > T > fcj, k we find 



-in 
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Im n T (p , p) = -— / dk / du[(k-w)/3 + (w,k) + (k + w)0-(v,k)], (6) 



-ft 



where [3± = —ImD± /it are the discontinuous parts of the spectral functions belonging to the effective fermion 
propagator. This expression agrees up to some prefactors with the photon production rate. Thus following Kapusta 
et al. |33|] we end up with 

The hard contribution, for which the momentum of the exchanged fermion is larger than A, can be calculated most 
conveniently from the matrix elements according to Fig.l of Compton scattering 

iZTI = j~J n F {k) J [1 + n B (p')] J [1 - n F (k>)] 

(2tt) 4 «5 4 (P + K-P'-K')4 (\M\ 2 ) comp (8) 



and pair creation 



TEH = Tp J ns(fc) / [1 nF[P ' )] J (2^F [1 - nF(k ' )] 

(2^) 4 S\P + K-P'-K')2 (\M\ 2 ) pair . (9) 
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Here momenta with a prime belong to outgoing particles. The factors in front of the amplitudes, averaged over initial 
states and summed over final ones, come from summing over the possible states of the incoming thermal particle. The 



matrix elements using Mandelstam variables, s = (P + K) , t = (P — P') 



'\2 



-s — t, are given by 



(\M\ 2 ) comp ^-2e i (- + - 



(\M\ 2 ) pair = 2, 



(10) 



In order to evaluate the integrations over the final states in (||) and (||) we assume p' T and k' ^> T, which 

is justified since p > p' + k! ^> T and the phase space for p'^T or k'^T is unfavorable j53|. Thus we may assume 
1 + 71b{p') — l~np(p') ~ 1 — np(k') ~ 1 corresponding to the Boltzmann approximation for the distribution functions, 
e.g. 1 + ub(p') — exp(p'/T) nsip') — 1, also used by Kapusta et al. |Q in the case of photon production, where, 
however, an integration over the initial states was considered. This approximation simplifies the expressions (^|) and 
(^J) considerably because we can now evaluate the integrals over p' and k' transforming to the center of mass system 
using the Lorentz invariant phase space factor |3~|| 



dQ = (2tt) 4 6 4 (P + K-P'- K 1 ) 



d 3 p' 



d 3 k' 



dt 



(2tt) 3 2p' (27r) 3 2fc' 8tts' 



Let us first consider Compton scattering: 



J dQ(\M\ 2 ) comp = ^ ( ln ^ + ^)' 



(11) 



(12) 



where A 2 <C s cuts off the logarithmic infrared divergence of the u-channel. 
Using s — 2pk(l — p • k) we obtain after integrating over d 3 k 



comp 
^hard 



_7r_ a 2 T 2 
12 p 



In 



8pT 



■7 + 



TO 

C(2) 



(13) 



where 7 = 0.57722 is Euler's constant and Q{z) is Riemann's zeta function with C / (2)/C(2) = —0.56996 
The pair creation contribution gives analogously 



pair 
^hard 



7T a 2 T 2 
6 p 



In 



4pT 
~A2~ 



1-7 + 



TO 

C(2) 



Adding up the ([l3|) and ( |T^ ) we get the total hard contribution 



1 a2T2 
4 p 



In 



pT 1 

A2 + 3 



In 128 



7 + 



TO 

C(2) 



(14) 



(15) 



Adding up the soft contribution (Q) and the hard one (|T^) the separation scale A cancels, as expected, yielding the 
final result 



7= 4 



7T a 2 T 2 , 0.3090p 



In ■ 



aT 



(16) 



Inserting into this expression a typical temperature for a supernova core, T — 10 MeV, (neglecting a finite chemical 
potential) a 100 MeV 7-ray has a mean free path A = 7~ 1 /2 of 0.40 nm. 

In the case of a quark-gluon plasma the photon damping is caused by the processes of Fig.l, where the outgoing 
photon for Compton scattering and one incoming for pair creation has to be replaced by a gluon, leading to 



7 = T 



5tt aa s T 2 , 0.2317p 



In- 



P 



(17) 



The mean free path of a 1 GeV photon in a QGP at T = 200 MeV, extrapolating the result ( p~7| ) obtained in the 
weak coupling limit g <C 1 to a realistic a s = 0.3, is A = 480 fm. This value is much larger than the dimensions of a 
fireball in ultrarelativistic heavy ion collisions confirming that photons may be used as a direct probe of the fireball 
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J33| , p5[] . (Of course, higher order processes such as bremsstrahlung may contribute to the rate similarly for realistic 
values of the strong coupling constant.) 

Comparing (17) with the photon production rate pdR/d 3 p |33|,|36| we observe that the latter differs only by a 
factor 4/(27r) 3 exp(— p/T). The Boltzmann factor reflects the principle of detailed balance relating the decay rate to 
the inverse rate 1^ by Td(p) — exp(— p/T) Ti(p) (3^]. The factor 4/(27r) 3 is just a matter of definition: pdR/d 3 p — 
2r.;/(27r) 3 and 7 = 1^/2. This also shows that the hard part of the photon production rate |]33| , j37f can be calculated 
much easier by considering the decay instead of the production process. 
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FIG. 1. Lowest order Feynman diagrams for photon damping: (a) Compton scattering, (b) pair creation. 



FIG. 2. Photon self energy diagrams containing an effective fermion propagator denoted by a blob. 
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This figure "figl-l.png" is available in "png" format from: 
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